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It has been shown that well-behaved spacetimes may induce the vacuum fluctuations of some 
nonminimally coupled free scalar fields to go through a phase of exponential growth. Here, we 
discuss this mechanism in the context of spheroidal thin shells emphasizing the consequences of 
deviations from spherical symmetry. 
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I. INTRODUCTION 



In a recent paper, it was shown that certain well- 
behaved spacetimes are able to induce an exponential 
enhancement of the vactium fluctuations of some non- 
minimally coupled free scalar fields [l|. This "vacuum 
awakening mechanism" may have consequences, in par- 
ticular, to astrophysics, since the vacuum energy density 
of the scalar field can grow as large as the nuclear density 
of neutron stars in few milliseconds once the effect is trig- 
gered [^]. As a result, the system must evolve into a new 
equilibrium configuration and eventually it should induce 
a burst of free scalar particles Q (see also Ref. [J|-[3 for 
related classical analyses reaching similar conclusions). 
Conversely, the existence of classes of nonminimally cou- 
pled scalar fields can be unfavored by the determination 
of the mass-radius ratio of relativistic stars with known 
equations of state. 

It is, thus, interesting to know if the main features 
described in Ref. are preserved when assumptions as 
staticity and spherical symmetry are relaxed. In this pa- 
per, we investigate the vacuum awakening mechanism in 
the context of thin static spheroidal shells. This will al- 
low us to explore the consequences of deviations from 
sphericity, while avoiding complications concerning un- 
certainties about the interior spacetime of nonspherical 
compact sources. 

The paper is organized as follows. In Sec. |lll we fol- 
low Ref. [6] and present the general properties of the 
shell spacetime, emphasizing the assumptions which were 
made in order to obtain the particular class of solutions 
that we investigate. In Scc. lIIIl we consider the quantiza- 
tion of a real scalar field in this background and proceed 
to discuss the vacuum awakening effect in nonspherical 
configurations. We show, in particular, that in the limit 
where spherical symmetry is recovered our results can be 
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expressed in terms of known functions. In Sec. lIV| we dis- 
cuss the exponential growth of the vacuum energy density 
in the context of spherically-symmetric shells. Section IVl 
is dedicated to conclusions. We assume natural units 
in which c = h = G - 1 and metric signature (- + ++) 
throughout the paper. 



II. THIN SPHEROIDAL SHELLS 

Let us consider a static and axially-symmetric thin 
shell surrounded by vacuum |7|-[8|. The most general 
line element describing the external- and internal-to-the- 
shell portions of the spacetime cornplying with the as- 
sumptions above can be written as [9| 

ds^ = -e^^dt^ + e2(--^)(dp2 + ^^2) ^ p^e-^^dip^, (1) 
where A = A(p, z) and v = v{p, z) satisfy 
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dp 
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(4) 



The external- and internal-to-the-shell regions will be 
covered with coordinates {t,p,z,ip) and {t,p,z,Cp), re- 
spectively, where we will denote by iS the 3-dimensional 
timelike boundary between them. It is worth to note 
that by using the spacetime symmetries one can choose 
the time and angular coordinates on S such that i = t 
and (p = Lp. As a result, we will denote the internal co- 
ordinates simply as (i, p, z, (p) and the shell is identified 
with t = const sections of S. Let us assume, moreover, 
that the shell lies on a A = const surface: 



X{p,z)\s = \{p,z)\s = Ao = const. 



(5) 



This choice leads the spacetime inside the shell to be flat 
with the corresponding line element being cast as jfij 

dsl = -e^^°df + e'^^" {dp^ + dz^ + p^d^^). (6) 

In order to analyze the external metric, it is convenient 
to perform the coordinate transformation {p, z} ->■ {x, y} 
defined by 



pEa(a;2- 1)1/2(1 -y2)i/2^ ^^ 



axy, 



(7) 



where x € [l,oo), y e [-1,1], and a = const > 0. In terms 
of the X and y coordinates, Eq. ([2]) reads 



d_ 
dx 



ox 



d_ 
dy 



2.dX 



dy 



0. 



(8) 



The most general solution of Eq. ([8]) which is regular 
on 2/ = ±1 (symmetry axis) and well behaved at x ^ oo 
(spatial infinity) can be cast as 



A = ^ AjQj{x)Pj{y), Aj = const. 



(9) 



3=0 



where Pj{z) and Qj{z) are the zero-order associated Leg- 
endre functions of first and second kinds [l3], respec- 
tively. For the sake of simplicity, we will restrict ourselves 
to the particular class of spheroidal shells obtained by im- 
posing Aj = for j = 1,2,... in Eq. ©. As a result, we 
have 



A: 



/3, x+1 

--log -, 

2 X - 1 



(10) 



where /? = -Aq > will play the role of a geometric pa- 
rameter linked to the shell shape. It is worthwhile to note 
that condition ([5]) combined with the Ai = A2 = . . . = 
choice restrict the possible shapes and stress-energy- 
momentum distributions of the shells considered here 
(see, e.g., Ref. [11] for more general shells). Still, this 
class of shells is general enough for our purposes. 

Equation (|T0)) implies that the shells which we consider 
will lie on 

X = xq = const > 1 

surfaces. The corresponding v solution can be directly 
obtained from Eqs. ©-Q: 



f3\ x2-l 
v = — log- 



(11) 



2 "°a;2-y2- 
By combining these results, the exterior metric will read 

*a^[^)'(x'-l){l-p')df,'. (12) 



One can see that the spacetime is asymptotically flat by 
taking the a; ^ +00 limit in Eq. (1121) . It will be shown 
later that 0</3< 1, (3 = 1, and /3 > 1 are associated 
with prolate, spherical, and oblate configurations, respec- 
tively. 

Next, we must impose continuity of the internal and 
external induced metrics, hab, on 5. It is convenient to 
cover S with coordinates ("" = (t^y^'P), « = 0,2,3, since 
the shell lies at a; = xq = const. The continuity condition 
establishes a relationship between the internal, p, z, and 
external, y, coordinates on S. For further convenience, 
however, let us replace coordinates p, z by f , as defined 
below: 



ar sm f 



ar cos ( 



By doing so, Eq. ([6]) reads 



ds^--i^-^'de.aH^-^' 
Vxo + l/ Vxo-l/ 

X (df2 + f2^Q2 ^ ~2 gj^2 g^^2-^ ^ 



(13) 



where we have used Aq = -(/3/2) log[(a;o + l)/(a;o - 1)] 
Then, in order to join the metrics given by Eqs. (fT2] 
and (|13|) on S, we impose f\s = f{y) and cos^jg = g{y) 
where 



and 



f{y) 



9'{y){i-y') 
i-9iyr 



{xl-l){l- y') 

^-giyy 



.1/2 



(14) 



5(y)2/+(l-.9(2/) ) 



2U/2 



\^l-y^) 



p^-i 



-.1/2 



(15) 



with " ' " E d/dy. The f{y) function follows immediately 
after g{y) is determined from Eq. (jlSp through a numeri- 
cal calculation, where we fix g{0) = (to harmonize with 
the solution g{y) = y for /3 = 1). In order to guaran- 
tee that g{y) is real, an extra restriction on xq must be 
imposed when /? > 1: 



/3 > 1 ^ xo > /?, 
0<^<l^xo>l. 



(16) 
(17) 



For /? = 1, Eq. (ITTI) just reflects the fact that the radius 
of a spherical shell must be larger than the Schwarzschild 
one. 

In order to investigate the dependence of the shell 
shape on /? and xq, we calculate 



^ — -^equatorial/ -'^meridional 
(l_l/a;2)(l-/3^)/2 

%i^i(l/2,(/32-l)/2,l;Xo2)' 



(18) 
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FIG. 1: The ratio n = I/equatorial /^meridional IS plotted as a 
function of /3 for shells lying at different values of xq = const. 
For < /3 < 1, /3 = 1, and /3 > 1 we have prolate (0 < n < 1), 
spherical (n = 1), and oblate (n > 1) shells, respectively. The 
maximum value of n, Wmax ~ 1.3, corresponds to an oblate 
configuration lying at a; = xo = /? ^ oo for which the equatorial 
diameter is about twice as large as the polar one. Shells with 
n ^ correspond to infinitely thin and long rods, < /? < 1, 
lying at a; = a;o -* 1. 



where ^meridional and ^equatorial are the meridional {ip = 
const) and equatorial (y = 0) proper lengths, respectively, 
taken on some t = const hypersurface and 2Fi{a,b,c; z) 
is the hypergeometric function. We note that the shell 
will be prolate (0 < n < 1), spherical (n = 1), and oblate 
(n > 1) for < /3 < 1, /3 = 1, and /3 > 1, respectively 
(see Fig. [T|) . The maximum oblateness (associated with 
the maximum n value) can be obtained from Eq. (11^ 
combined with Eq. (IK 



lim 



=.1/4 



//o(l/4) « 1.3, 



where Iiy{z) is the modified Bessel function of first kind. 
This limit should not be viewed as a general restriction 
to arbitrary oblate shells but rather a consequence of 
the assumptions discussed below Eq. pI7)) . There is no 
similar restriction for prolate configurations, since n may 
take arbitrarily small (positive) values. For the sake of 
further convenience, it is also useful to calculate, at this 
point, the shell proper area as a function of (3 and a;o: 



A- 
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(;3" + l)/2 



(20) 



By establishing the interior and exterior metrics, the 
shell stress-energy-momentum tensor is also fixed: 



j.p.^^a6gMg^5(^)^ 



(21) 



where £ is the proper distance along geodesies intercept- 
ing orthogonally S (such that £ < 0, i = 0, and f > in- 
side, on, and outside S, respectively), e^ = dx*^ jdC,"" are 



the components of the coordinate vectors d/d^"" defined 
on S (with {x^} being some smooth coordinate system 
covering a neighborhood of 5 [S]), and 



gab 



8n 



h''''AK). 



(22) 



Here, Kab is the extrinsic curvature, K = Kabh"^^, and 
AA"*^„Q gives the discontinuity of ^°m„o... across S. 
A straightforward calculation leads to f(|] 



87r5"o = A{y)[B{y) + C{y) - 2p{xl - y^)], 
STrSl = A{y)C{y), 



87TSl = A{y)B{y), 



(23) 
(24) 
(25) 



where 



A{y) . ^, 
a{x, 



-l)-^lxo-lf'(xl-y'\ 
l-y2)\^, + l) [xl-lj 



,(r-i)/2 



B{y) E l3'xo{l - y') - U{y)-\xl - 13Y){4 " 1) 
+ Xo{xl-l), 

Ciy)^ixl-y')[xo-Uiy)l 



and 



U{y) 



\l-y^) [xl-y^J 



V 



1/2 



(26) 



Then, by using Eq. ^ combined with Eqs. (I23 l) -p5 l) . 
we obtain that the gravitational mass formula [12] 



M = 2 



ih-h'^^'y'''" 



(19) yields 



M = Ba. 



(27) 



(28) 



Here, <j'' = {d/dtY is a global timelike Killing field, the 
integral is taken on a i = const Cauchy surface Et, and 
dE^ = Ufj^dH with n'^ being the pointing-to-the-future 
unit vector field orthogonal to Et. 

It can be verified that the weak and strong energy 
conditions are always satisfied by the stress-energy- 
momentum tensor (j2ip . As for the dominant-energy con- 
dition, it will be satisfied for < /3 < 1 and /3 > 1 provided 
that 

2^x0 > -1+13^ + 3x1- 4' (xq - 1)^"^'/^ 

-2xl^\xl-ir'^ (29) 

and 

xo>(13/12)/3, (30) 

respectively (see Fig. [5]). (One can see that for /3 = 1 
Eqs. (l29l) and (|30)) agree with each other.) Therefore, 
the matter composing this class of shells has reasonable 
physical properties for a significant range of parameters. 




FIG. 2: The light gray region corresponds to values of {/3,xo) 
which do not satisfy the dominant-energy condition, while the 
dark gray region is excluded by the constraint (|16|) . The blank 
area corresponds to shell configurations satisfying the weak, 
strong and dominant energy conditions. 



III. QUANTIZING THE FIELD AND AWAKING 
THE VACUUM 

Now, let us consider a nonminimally coupled real scalar 
field $ with null mass defined over a spacetime of a 
spheroidal shell as discussed in Sec. |ll] (see Sec. Ill of 
Ref. 3] for a discussion about the physical reasonable- 
ness of the null-mass assumption). It will satisfy the 
Klein-Gordon equation 



VuV''$ + C^$ 



0, 



(31) 



where R is the scalar curvature and ^ = const is a dimen- 
sionless parameter. 

We follow the canonical quantization procedure and 
expand the corresponding field operator as usually: 



$: 



d'd{r])[ajjUjj + a^u^], (32) 

where i9 is a measure defined on the set of quantum num- 
bers rj. Here, m^ and u* are positive and negative norm 
modes with respect to the Klein-Gordon inner prod- 
uct [13], respectively, satisfying Eq. (PT|) . Then, the an- 
nihilation a^ and creation aj operators satisfy the usual 

commutation relations [aQ,,aI] = 5{a,f3), [0^,0^] = 
and the vacuum state |0) is defined by requiring a^|0) = 
for all rj. 

Since the spacetime is static and axially symmetric, it 
is natural to look for positive-norm modes in the form 



,(i,x,^) = r,(0^^p(x)e^'^^ 



(33) 



where x - (^:^) ^^'^ X - (^jV) inside and outside the 
shell, respectively, /i 6 Z is the azimuthal quantum num- 
ber, and a = const. By using Eq. ([55)1 in Eq. (PT)) . we see 
that Trj{t) obeys 



d^ 



0, 



(34) 



while FcTf^ix) satisfies 



1 rxo-l \ 



2,3 



,2 






deism 9 de) 



f^ sin ( 



K 



<Tfl 



<jF. 



(TfJ, 



(35) 



and 



a2 U + 1 / 



2/3 



ix'-y') 



■2\0^-l 



(a;2-l)^ 



■{d4{x^-l)d,] 



dy[i^-y')dy]}- 



(x2-l)(l-y2) 



F, 



(Tfl 



<^Kf.- 



(36) 



Here, we have assigned labels "-" and "+" to F^.^ in 
order to denote solutions valid inside and outside the 
shell, respectively. 

The solutions of Eq. (|34p will assume the following gen- 
eral forms: 



T„it)^T^it)ocexp{-iujt), w>0, 
for a = Ld^ > and 



+ e 



(37) 



n>0, (38) 



for a = -il^ < 0, where the latter is one of the possi- 
ble combinations which guarantee that u„ with cr < 
is indeed a positive- norm mode [i|. Equation ([37]) is 
connected with the usual time-oscillating modes while 
Eq. ([38t is associated with the so-called "tachyonic" 
modes. Tachyonic modes are responsible for an exponen- 
tial growth of quantum fiuctuations and, consequently, 
of the expectation value of the stress-energy-momentum 
tensor [V\. (We address to Ref. [3| for more details on 
the canonical quantization procedure in the presence of 
unstable modes but it is worthwhile to emphasize at this 
point that tachyonic modes do not violate any causality 
canon.) The requirement that these modes be normaliz- 
able determines the possible negative values of a (if any) 
and, thus, the existence (or nonexistence) of tachyonic 
modes. One sees from Eq. ([55)1 that tachyonic modes 
vanish exponentially at infinity: 



K^i^iV) -* Fn^i^iV) ''S:°°exp(-r2c 



(39) 



Next, let us analyze Eqs. ((35)) - (|36l) in more detail. On 
account of Eq. (PTj) . we have that 

i?=-87rT=-2AX(5(£). 

Then, one sees from Eq. (pTj) that F^ should join each 
other continuously on S: 



ail 



(f,0)|5=F;^(x,y)|< 



(40) 



while the first derivative of -F^^ along the direction or- 
thogonal to the shell will be discontinuous: 



^{dF„iild(.)\s = ^"t{v)F,^\s. 



(41) 



Here, j{y) = -2AK and we recall that AK = 47r(S'°o + 
5 2 + S%). By using Eqs. ([231) -([25]), we obtain 



2 ( XQ-I Y ( xl-y^ y 



P'^Xo - P + Xq 



^0-1 



, xo(i-/3^) 1 / yHP'-i) \ U{y) 

" " Uiy)\ xl-y^ ) xl-1 



a;2-y2 



(42) 



where we recall that U{y) is given by Eq. I^E^. For the 
sake of convenience, we cast Eq. (1411) in a more explicit 
form: 



\di dx j g [ 



di dx I g \d(. df ^ d£ 89 



(43) 



where 



dx 

H 



s a\xQ + ll \xl-l ) 



2\(r-i)/2 



and 



with 



df 

de 



d£ 



a \xo + l / 



lx/3/2 



D{y)9'{y) 
\/i-5(y)^' 



D{y)ny) 
fiyy 



9'{y? . fiyy 



1-1/2 



i-9{yr /(y)M 



(44) 



(45) 



(46) 



(47) 



In what follows, we search for the ^ parameters which 
give rise to tachyonic modes and, hence, to the vacuum 
awakening effect, once the spacetime, characterized by 
the values of xo, (3 and M, is fixed. For this purpose, we 
must look for regular F^ (x) functions with ct < satisfy- 
ing Eqs. ([55]) and ([55]) inside and outside the shell, respec- 
tively, while respecting Eqs. (|40|) and (|4T|) on the shell 
and vanishing exponentially at infinity [see Eq. (139^ ] . 



A. Spherical shells 

Let us start with an analytical investigation of the con- 
ditions required by spherically-symmetric shells to allow 
the existence of tachyonic modes. It will be interesting in 
its own right and useful as a test for the reliability of the 
numerical code which will be used to treat more general 
axially-symmetric shells further. 

First, we note from Eq. ^^ that a = M for /3 = 1. 
Hence, by using the definitions x = r/M - 1 and y = cos6 
in Eq. ^ and xq = R/M - 1 in Eq. ^, we write the 
external- and internal-to-the-shell line elements as 



and 

ds^_ = -(1 - 2MIR)di^ + (1 - 2M/R)-iM2(df2 + f^ds^), 

(49) 
respectively, where ds^g = dO^ + sin^ Odip^ . In terms of the 
internal and external coordinates, the shell is at 

f = (1 - 2A///R)^/2R/Af and r = R, 

respectively, from which we see that R is indeed the shell 
proper radius. In order to define a continuous radial co- 
ordinate on the shell, we introduce 

r_EMf/(l-2M/R)^/^ r+ = r 

with respect to which the shell will be at r± = R. 

Now, we note that the general solutions of Eqs. ((35)) - 
([55)) can be cast in the form 

F^^{r^,0) = f a,^(V'^,(r,)/r,)P;(cos0), (50) 

where a;^ = const and Pj^{cosO) are associated Legendre 
functions of the first kind, degree Z = 0, 1, 2 . . ., and order 
/i= -Z,-Z + l,...,/. For the sake of convenience, we define 
the coordinates 

X_ = r_l{l-2MlRfl\ x+ = r^ + 2M\n[r^l{2M)-l] + D, 

where D = const is chosen such that x- ^^^id x+ fit each 
other continuously on the shell. By using x±i th^ func- 
tions ijj^i will satisfy the "Schrodinger-like" equation 



■ d^i^tildxl + Kff 'V^/ = 'J^ti 



with 



K 



(i,-) 



off 



(l-2M/R)/(/+l)/r2, 



yj^'^) = (1 - 2M/r,mi + l)/rl + 2M/rl). 



(51) 

(52) 
(53) 



The discontinuity of the potential across the shell is 
AV^ff = 2M(1 - 2M/R)/R3. Ahhough v}g^'^ is positive 
everywhere off shell, the existence of tachyonic modes is 
still possible because the potential on the shell contains 
a delta distribution. Thus, depending on R/M, ^ and /, 
the effective potential will be "negative enough" to allow 
solutions of Eq. (|5T)) for cr < 0. This is codified in the 
discontinuity condition (1411) . which will be used later. 

As a result of Eq. (|5(I)) . the field-operator expan- 
sion (l32l) can be written in this case as 



J, (54) 



Ifl 






nifj.''^nitj.h 



dsl = -{l-2M/r)dt^ + {l 



2M/ry^dr'^ + r'^dsl 



(48) 



where the only nonzero commutation relations between 
the creation and annihilation operators are 

[bujit„bl,i,^,] = Sw5f,f,,S{uj-u}'), (55) 

[cnii,, cli'i'n'^ = SwSf,f,'dnn', (56) 



the modes read 



vti^ = T^{t){r^i{T,)lr^)Yi^{9,v), 



WQin = Tn{t){ipni{r±)/r±)Yif,{e,ip) 



(57) 



(58) 



with 



Yi^ie,^) 



{2i + i){i-fiy. 

\| 4:7r{l + fi)l 



Pj" {cos e)e 



ifiip 



and we have assigned labels "±" to w^^;^ and w^ifi to de- 
note solutions valid inside and outside the shell, following 
our previous notation. 

Our search for tachyonic modes equals, thus, the search 
for solutions Vn; of Eq. ^^ with a = -fl'^ < 0. The 
regularity requirement for the normal modes implies that 
at the origin 






0\ 



(59) 



where we have assumed that it approaches zero from pos- 
itive values, since the Klein-Gordon inner product fixes 
the mode normalization up to an arbitrary multiplica- 
tive phase (which can be chosen at our convenience). By 
using Eq. (IFTj) with Eq. (15^ . we conclude that 



0<i'm{r-)\s = ^ni{r.)\i 



where the equality is a consequence of Eq. (^(11) . 

On the other hand, Eq. (|41l) implies that the first 
derivative of the radial function will be discontinuous on 
S: 



R / dr+ dr^ 



I 2M 



4 



"'\ ^niir-) 



(61) 



Then, dip^i/dx+ls will be, in general, a nontrivial func- 
tion of the field and shell parameters. Now, by noting 
from Eq. ([ST]) that 

^^ilO^ d^i^tJdxl I 0, 

we conclude that either tjj'^i changes sign once, diverg- 
ing negatively at infinity, or it remains always positive. 
Tachyonic modes (with il > 0) will be associated with 
TpQi > with the additional requirement that 



lim ijtii = 0^- 



(62) 



It follows, then, from Eq. (ISTI) that for a given shell con- 
figuration there will exist up to one tachyonic mode for 
each fixed I. 



In order to investigate which shell configurations give 
rise to tachyonic modes, let us first note that there always 
exist a negative enough a = -fig < 0, such that 



lim ipLi = +°°- 



X+^ + oo 



Then, if 



lim tp^i = -oo. 



X+^ + oo 



(63) 



(64) 



it is certain that there will exist some fl e (0,51o) satis- 
fying condition (|5^ . Conversely, if condition (|M| is not 
verified, there will be no tachyonic mode. 

The solutions of Eq. JST]) with ^ = satisfying Eq. ([59]) 
can be written as 



%i{r-) = AirJ^'/R\ Ai > 



(65) 



V'^,(r+) = BiPiirjM - l)r+ + CiQi{rjM - l)r+, (66) 

where Ai,Bi and Ci are constants. Next, by imposing 
conditions (pOj) and (pT|) on the shell, we obtain 



A, 
and 



- = !-— 2 + 2a/i In] 

m\r V R / V 



2M 



R ) ("' 



(60) Bi^ 



Bi _ [(/ + 40(Af.To/R-Vl-2M/R)-2M$/R]Q;(xo) 
Ai (IM/R) [Pi{xo)Qi-i{xo) - Pi-i{xo)Qi{xo)] 



{lM/R)Qi^i{xo) 



(IM/R) [Pi{xo)Qi-i{xo) - Pi-iixo)Qi{xo)] 

for / = and ^ > 1, respectively, and 

Q l-iBi/Ai)Piixo) 



Ai 



Qi{xo) 



(68) 



(69) 



for / > 0, where we recall that in the spherical case xq 
R/M - 1. Then, by using that 



lim Qi{rjM-l)r^ 



lim Pi{r+lM-l)i 






(70) 
(71) 



we see from Eq. (p5|) and condition (|M)) that the existence 
of tachyonic modes with some J7 > requires Bi/Ai < 0. 
{Bi/Ai = corresponds to "marginal" tachyonic modes 
characterized by having the quantum number Q = 0.) 
This establishes a relationship between the field parame- 
ter ^ and the shell ratio M/R. In Figs. [3] and SI we show 
the parameter-space region where tachyonic modes with 
/ = and ^ = 1 do exist, respectively. We note that be- 
cause the smaller the I the lower the V^fj'"^ , the existence 
of a tachyonic mode with I = Iq implies the existence of 
tachyonic modes with 1 = 0,...,Iq-1. This can be seen in 
Figs. [3] and 0] as we note that the tachyonic- mode region 






FIG. 3: The black and dark-gray areas depict the parameter- 
space region where Bo/Aq < leading to the "vacuum awak- 
ening effect" . The magnitude of the vacuum energy density 
on the shell grows positively and negatively in the black and 
dark-gray regions, respectively (see discussion in Sec lIVf) . The 
light-gray strip is excluded from the parameter space because 
no static spherical shell can exist with R < 2M, while the 
translucent-gray one contains those configurations which vi- 
olate the dominant-energy condition. The inset graph em- 
phasizes that there are shell configurations which allow the 
presence of tachyonic modes for ^ = 1/6 (vertical dashed line), 
although the dominant-energy condition is not satisfied. 



FIG. 4: The black areas depict the parameter-space region 
where tachyonic modes with I = 1 are present. The light- and 
translucent-gray regions represent the same as in Fig. |3] In 
contrast to the I = case, no analysis is performed here con- 
cerning whether the vacuum energy density on the shell grows 
positively or negatively because any contribution coming from 
/ = 1 is dominated by the one associated with I = 0. 



modes are 



Foo = ^Azf'Pz(cos( 



(72) 



for Z = 1 is contained in the one for 1 = 0. Clearly, the 
existence of a single tachyonic mode is enough to induce 
an exponential growth of quantum fluctuations leading 
to the vacuum awakening effect. We note, in particular, 
that there are shell configurations which allow the ex- 
istence of tachyonic modes for the conformal field case, 
f = 1/6. Nevertheless, it can be also seen from Figs. [3] 
and m that for these configurations the dominant-energy 
condition ([50]) is violated. 



B. Prolate and oblate shells 

Now, we proceed to treat the prolate (0 < /? < 1) and 
oblate {(3 > 1) spheroidal shell cases. Here, we shall focus 
our attention on the boundaries which curb the regions 
where the vacuum awakening effect occurs due to the 
existence of any tachyonic mode. These boundaries are 
associated with the presence of marginal tachyonic solu- 
tions with i7 = [see discussion below Eq. ([7T|)]. More- 
over, following the spherically-symmetric case reasoning 
where the most likely tachyonic modes have I = (im- 
plying /i = 0), we will look for marginal tachyonic modes 
(fl = => (7 = 0) with /i = in the axially-symmetric 
prolate and oblate cases. Then, the relevant regular solu- 
tions of Eqs. ([55)1 and ([55]) which give rise to normalizable 



and 



F^o = TCiQi{x)Pi{y). 



(73) 



We note that Eqs. ([7^ - (|75|) generalize the spherically- 
symmetric relation (j50l) with cr = ^ = provided that one 
sets Bi = in Eq. (p5)) [see observation within parentheses 
below Eq. ([7T1) ]. 

Next, we use Eqs. ((72|) - ((73)) in the continuity condi- 
tion (|1U)) to determine the C; coefficients in terms of the 
Ai ones: 

?/' + 1 ri 

where we recall that /(y) and g{y) are given in Eqs. (|14p 
and (jl5|) . respectively, and we have used the orthonor- 
mality condition 



X 



1 2 

Pn{y)Pm{y)dy= -Snm- (75) 

1 2n+ 1 



The boundaries which curb the regions possessing 
tachyonic modes are obtained as we impose the first- 
derivative constraint (j4ip . In order to obtain them, 
it is convenient to note that Eq. (pij) can be cast as 



Y,i AiGi{y) = for an intricate but otherwise known func- 
tion Gi{y). By expanding Gi{y) in terms of Legendre 
polynomials, Eq. (PT|) can be written as 



^^A,fc„,PKy) = o, 



(76) 



; ;' 



where kw = {I' + l/'2){klp + kj^ + kl^) with 






l(2) 



1 dQi,{x) 
Qi'{x) dx 



■ j\Pv{y) 

ifivY 



dx 



T-l 



f\Pv{y)f{yyPi[9{y)l 

(77) 
dr' 



if(y) 



i-i 



Pildiy)] 



VT^ffCyP 



d9 

de 



{-g{y)Pi[g{y)] + Pi-i[9{y)]} 



and 



(3) 






-ij\yPv{y) 



dx 

'de 



(78) 



i{y)f{yyPi[9{y)l (79) 



Here, we recall that 7(2/), dx/de\s, dr/di\s, and d6/de\s 
are given in Eqs. (H^ . (Hil) . PSI) . and (|i51) . respectively. 
Then, by using the orthonormality property of the Leg- 
endre polynomials, Eq. ((75)) leads to 



J^Aiku 



0. 



(80) 



Let us, now, consider kw as elements of a matrix /C. In 
the spherically-symmetric case, /C is diagonal: kw = giSw 
with gi being constants depending on the shell, M/R, 
and field, £,,1, parameters. The borderlines associated 
with the regions containing tachyonic modes for each / 
are obtained by solving gi = for ,f as a function of M/R. 
In the absence of spherical symmetry, the corresponding 
borderlines inside which tachyonic modes exist can be 
obtained similarly by vanishing the eigenvalues of /C. The 
vanishing-eigenvalue condition can be imposed on /C by 
solving the corresponding characteristic equation 



det/C = 0. 



(81) 



This will drive Eq. (|5(7)) to have a nontrivial solution for 
the Ai coefficients. We recall that, eventually, all modes 
should be Klein-Gordon orthonormalized, which fixes any 
remaining Ai left free. 

For computational purposes, we truncate (the infinite 
matrix) /C by imposing < ^,r < iV for a large enough 
N. This is justified since the ku/ elements decrease as 
the values of Z or T increase. By fixing the (3 and xq pa- 
rameters, Eq. (|5T|) is expected to be satisfied by A^ + 1 
values of £, which corresponds in the spherical case to 
the fact that for a fixed M/R value, the boundary of the 
unstable regions, associated with the marginal tachyonic 
modes, are at different ^ values; each one corresponding 
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FIG. 5: Diagram showing how higher order A'' approximations 
converge to the actual boundaries which curb the tachyonic 
unstable regions for a considerably prolate shell, n = 0.25. 
Configurations allowing for tachyonic modes are those to the 
left of the curves on the left-hand side and to the right of 
the curves on the right-hand side. We note that in the 
16tvM'^/A « 1 regime, the N = approximation is already 
quite satisfactory in contrast to the IGvrM /A » 1 regime. 
For N = 8, only few points were obtained due to the compu- 
tational cost. 



to a distinct I = 0,. . . ,N (see Figs. [3] and E}. Because 
in the prolate and oblate cases we are interested in the 
regions where the vacuum awakening effect occurs by the 
existence of any tachyonic mode, we shall look for the ^ 
solutions of Eq. (|5T|) which lead to the boundary enclos- 
ing the largest possible unstable region. (In the spher- 
ical case, it corresponds to the boundary of the black 
and dark-grey regions in Fig. [3] associated with I = 0.) 
For relatively small deviations from sphericity, a quite 
reasonable approximation is already obtained by taking 
N = 0. By increasing TV, we introduce higher order cor- 
rections. These corrections are seen to be more relevant 
as larger deviations from sphericity are considered as can 
be verified in Fig. [5] for a considerably prolate shell. 

In Figs. [6] and [71 we show the results obtained for some 
prolate and oblate shells, respectively, assuming iV = 6. 
For the sake of clarity, we have characterized the shells by 
their equatorial-per-meridional size ratios n and proper 
areas A as given in Eqs. (|18|) and ([20l) . respectively, since 
they have a more straightforward physical meaning than 
Xq and /?. Figure[6]shows that the lines which limit the re- 
gions where the vacuum is awakened by prolate shells dif- 
fer significantly from the spherical case for dense enough 
configurations, IdnM'^/A rj 1. In contrast to the spheri- 
cal case, where there is no equilibrium configuration for 
WttM'^/A > 1, in the prolate one, 16ttM^/A can acquire 
arbitrarily large values when n is arbitrarily small. Fig- 
ure [7] puts in context the oblate case. We recall from 
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FIG. 6: Diagram showing the boundaries which circumscribe 
the regions where the vacuum awakening effect is triggered 
by prolate-spheroidal shells with n = 0.25, 0.5 and 0.75. The 
spherical (n = 1) case is plotted for comparison. The vertical 
dashed line indicates the conformal-coupling value ^ = 1/6. 
Configurations allowing for tachyonic modes are those to the 
left of the curves on the left-hand side and to the right of the 
curves on the right-hand side. 



Sec. |IT] that the degree of nonsphericity for this class of 
solutions is restricted on account of the constraint (IT51) : 
/3 > 1 => xo > /3, leading to 1 < n < 1.3. This restriction 
reflects itself on the allowed values for WttM^/A. The 
excised regions at the top and bottom of Fig. [7] come 
from this condition applied to the oblate shell considered 
in the graph. From this figure, we also see that the oblate 
shell with n = 1/0.85 « 1.18 is more favorable to trigger 
the effect than the associated prolate one with n = 0.85. 



IV. EXPONENTIAL GROWTH OF THE 
VACUUM ENERGY DENSITY 

Finally, we investigate the exponential growth of the 
vacuum energy density induced by the existence of tachy- 
onic modes. Although the vacuum energy density will be 
a nontrivial point-dependent function, the total vacuum 
energy will be time conserved [l|. Let us suppose that 
a spheroidal shell evolves from (i) an initial static con- 
figuration {xo,l3,M)in where Eq. (PT|) is only allowed to 
have time-oscillating solutions to (ii) a new static config- 
uration (a;o,/3,M)out where Eq. (PT|) is permitted to also 
have tachyonic ones. Because the oscillating in-modes 
will eventually evolve into tachyonic and oscillating out- 
modes, the vacuum energy density (Too) - (0in|2oo|0in) 
will grow exponentially. Here, we assume the vacuum 
|Oin) to be the no-particle state defined according to the 
oscillating in-modes (see, e.g., Ref. 13| for a more com- 



FIG. 7: Diagram showing the boundaries which limit the 
regions where the vacuum awakening effect is triggered by 
an oblate-spheroidal shell with n = 1/0.85 ~ 1.18. Prolate, 
n = 0.85, and spherical, n = 1, cases are plotted, as well, for 
the sake of comparison (restricted to the domain of the oblate 
case, which is indicated by the horizontal dotted lines). The 
vertical dashed line indicates the conformal-coupling value 
5 = 1/6. Configurations allowing for tachyonic modes are 
those to the left of the curves on the left-hand side and to 
the right of the curves on the right-hand side. 



prehensive discussion). 

A general expression for the exponential growth of the 
expectation value of the stress-energy-momentum tensor 
was calculated in Ref. [ij . By applying it to the spherical 
shell case, we obtain the following leading contribution 
to the vacuum energy density: 



(Too) = (Too)H{-i) + (f^o)H{i) + (foo)s, 



(82) 



where we recall that £ is the proper distance along 
geodesies intercepting orthogonally S [as defined below 
Eq. (PT|) ] and H{£) is the Heaviside step function. Here, 
(Tqq), {Tqq), and (Too)5 are the vacuum contributions 
to the energy density inside, outside and on the shell, 
respectively, where 
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with K being a positive constant of order one related to 
the decomposition of the in-modes in terms of the out- 
modes and d, denoting the largest H, selected from the set 
of all tachyonic solutions. By analyzing the factor multi- 
plying the delta distribution in Eq. (j85p . one can verify 
whether the vacuum contribution to the energy density 
is positive or negative on the shell. Our conclusions are 
depicted in Fig. [3j the black and dark-grey regions are 
associated with shell configurations where {Too)s grows 
positively and negatively, respectively. Similarly, one 
concludes from Eq. (|83p that the total vacuum energy 
inside our spherical shells is positive, null and negative 
when ^ < 1/4, ^ = 1/4, and £_ > 1/4, respectively. 

As a matter of fact, eventually the scalar field and 
background spacetime must evolve into some final stable 
configuration, where tachyonic modes are not present, 
in order to detain the exponential growth of the stress- 
energy-momentum tensor. The precise dynamical de- 
scription of how the "vacuum falls asleep" again is 
presently under debate ^l^] . In spite of the quantum sub- 
tleties involved in this discussion, at some point the scalar 
field is expected to loose coherence after what a classi- 
cal general-relativistic analysis should be suitable. The 
evolution of axially-symmetric rather than spherically- 
symmetric systems may lead to new interesting features. 



V. CONCLUSIONS 

It was recently shown that relativistic stars are able to 
induce an exponential enhancement of the vacuum fluc- 



tuations for some nonminimally coupled free scalar fields. 
In Ref. [2| it was assumed spherical symmetry to describe 
compact objects, which is expected to be a very good 
approximation for most relativistic stars J15J. In this 
article, we were interested in analyzing how deviations 
from sphericity would impact on the vacuum awakening 
effect. For this purpose we have considered a class of 
axially-symmetric spheroidal shells. This has allowed us 
to pursue our goal, while avoiding concerns about how 
to model the interior spacetime of nonspherical compact 
sources. Fig. [5] shows that for dense enough configura- 
tions, 16ttM^/A Ri 1, the awakening of the vacuum be- 
comes more sensitive to prolate deviations from spheric- 
ity. Fig. [7] unveils that oblate shells with n = Uq seem to 
be more efficient to awake the vacuum in comparison to 
prolate ones with n = I/uq. 

As a consistency check, we have performed an ana- 
lytical investigation for the spherically-symmetric shell 
case in order to test the numerical codes used to dis- 
cuss the general axially-symmetric one. It was shown, in 
particular, that in contrast to the relativistic stars ana- 
lyzed in Ref. Q], spherically-symmetric shells are able to 
awake the vacuum for conformally coupled scalar fields, 
£_ = II & (see Figs. [3]|3]). The exponential growth of the 
vacuum energy density was analyzed for the spherically- 
symmetric case in Sec. IIVI 
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